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ABSTRACT 

The turbulent diffusion tensor describing the evolution of the mean concentration of a passive 
scalar is investigated for non-helically forced turbulence in the presence of rotation or a mag- 
netic field. With rotation, the Coriolis force causes a sideways deflection of the flux of mean 
concentration. Within the magnetohydrodynamics approximation there is no analogous effect 
from the magnetic field because the effects on the flow do not depend on the sign of the field. 
Both rotation and magnetic fields tend to suppress turbulent transport, but this suppression is 
weaker in the direction along the magnetic field. Turbulent transport along the rotation axis 
is not strongly affected by rotation, except on shorter length scales, i.e., when the scale of the 
variation of the mean field becomes comparable with the scale of the energy-carrying eddies. 
These results are discussed in the context of anisotropic convective energy transport in the 
Sun. 
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1 INTRODUCTION 

In the outer 30% of the sun, energy is transported mainly by con- 
vection through the exchange of fluid parcels. This process tends to 
smear out entropy gradients by making the gas nearly isentropic on 
the average. Any remaining entropy gradients lead to a mean en- 
tropy flux whose negative divergence leads to a local change in the 
entropy. This entropy flux is equal to the negative gradient of the 
mean entropy multiplied by some turbulent diffusivity coefficient. 
Using mixing length theory, this turbulent diffusion coefficient is 
found to be of the order of the turbulent velocity times some mix- 
ing length, which is often taken as a free parameter (Vitense 1953). 

With rotation, the smearing out of entropy gradients is sig- 
nificantly reduced due to lateral mixing generated by Coriolis de- 
flections (Julien et al. 1996). Such deflections cause fluid parcels 
to lose their heat content to the sides before they can ascend or 
descend and transport heat vertically. This additional mixing ex- 
plains why the smearing out of entropy gradients is heavily reduced 
by rotation. Rotational effects on the turbulent diffusivity were al- 
ready studied by Weiss (1965) and Durney & Roxburgh (1971) who 
found that the turbulent diffusion tensor becomes anisotropic with 
respect to the direction of the rotation axis. The components of 
the anisotropic diffusion tensor can be computed using the first- 
order smoothing approximation (Rudiger 1989, Kitchatinov et al. 
1994). These calculations show that the enhancement of turbulent 
diffusion along the direction of rotation corresponds to a latitudinal 
heat flux toward the pole, making it slightly hotter. Turbulent con- 
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vection simulations by Rudiger et al. (2005) confirm this, but they 
also show that the radial heat flux is slightly smaller near the poles, 
which is explained by the effect of stratificational anisotropy. This 
agrees with earlier results of Rieutord et al. (1994) and Kapyla et 
al. (2004). Rudiger et al. (2005) also established the presence of 
an azimuthal turbulent heat flux and found it to be negative, i.e., in 
the westward direction. This was surprising in view of the result of 
Kitchatinov et al. (1994) who found it to be either zero (for a sim- 
ple mixing length model) or positive (for a turbulence model with 
finite-time correlation effects). However, their result was based on a 
sign error (Rudiger, private communication). Furthermore, if there 
is also stratification, the simple mixing length model of Kitchati- 
nov et al. (1994) predicts that the sign of the flux should agree with 
that of the corresponding component of the Reynolds stress tensor 
and hence negative, which is then in agreement with the convection 
simulations of Rudiger et al. (2005). 

Diffusive processes also play a role in the evolution of the 
mean magnetic field and the mean momentum. In recent years it 
has become possible to compute the relevant diffusion coefficients 
together with other (e.g. non-diffusive) turbulent transport coeffi- 
cients using turbulence simulations. In the case of magnetic fields, 
a fairly accurate method is the test field method in which one solves 
an additional set of passive evolution equations for the departures 
from the corresponding set of prescribed mean fields (Schrinner et 
al. 2005, 2007; Brandenburg 2005, Brandenburg et al. 2008a). 

In cases of poor scale separation, i.e. when the scale of the 
mean field is only a few times larger than the typical scale of the 
turbulence, the multiplication with a diffusion coefficient must be 
replaced by a convolution with an integral kernel. A similar situa- 
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tion also applies to cases where the timescale of the variation of the 
mean fields becomes comparable with the turnover time. In such 
cases, the test field method has also been used to determine the in- 
tegral kernels (Brandenburg et al. 2008a, Hubbard & Brandenburg 
2008). If the turbulence is anisotropic and/or inhomogeneous, the 
diffusion coefficients or kernels are replaced by tensors. Again, the 
test field method allows an accurate determination of all relevant 
tensor components. 

In this paper, we apply the test field method to the turbulent 
diffusion of entropy in cases where the turbulence is homogeneous, 
but anisotropic due to the presence of rotation or a magnetic field. 
The former problem has recently been addressed by Riidiger et 
al. (2005), who evaluated the actual fluxes in stratified turbulent 
convection. However, convection necessarily imposes an additional 
anisotropy in the direction of gravity, so it becomes harder to disen- 
tangle the effects of gravity and rotation. Here, we focus on the case 
of forced turbulence in the absence of gravity; using an isothermal 
equation of state. In that case, the evolution equation of entropy is 
similar to that of a passive scalar. 

The effect of rotation on turbulent diffusion is important for 
modeling the possibility of latitudinal entropy gradients in rotating 
stars, i.e. perpendicular to the radial direction of the energy flux. 
Latitudinal entropy gradients could be important for producing a 
significant baroclinic term that would balance the curl of the Cori- 
olis force and could therefore be important for explaining the ob- 
servation that the contours of constant angular velocity in the Sun 
are not parallel to the rotation axis (Kitchatinov & Riidiger 1995; 
Thompson et al. 2003). The effect of a magnetic field on turbulent 
heat diffusion is important in the theory of sunspots (e.g. Kitchati- 
nov, & Mazur 2000) where the energy transport is suppressed in 
an anisotropic fashion (Riidiger & Hollerbach 2004). Furthermore, 
anisotropic passive scalar transport is important in connection with 
understanding the partial lithium depletion in the Sun (Riidiger & 
Pipin 2001). 



replaced by a convolution with an integral kernel. We discuss this 
case near the end of the paper, but ignore this for now. 

The turbulent diffusion tensor is a proper tensor, so in the pres- 
ence of rotation and in the absence of helicity it has only even pow- 
ers of the pseudo vector $7 and odd powers of f2 combined with odd 
powers of the Levi-Civita tensor. It must therefore have the general 
form (Kitchatinov et al. 1994) 



(4) 



where f2 = f2/|n| is the unit vector along the rotation axis and 
Kt, kq, and kqq are functions of the Coriolis, Peclet, and Schmidt 
numbers, that are defined by 
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respectively. Incidentally, we note that kq contributes an anti- 
symmetric component to the turbulent diffusion tensor, which does 
not contribute to divergence of flux. Similarly, in the presence of 
a magnetic field B with unit vector B — B /\B\, the diffusion 
tensor must have the form 



(6) 



but here kb = 0, because the diffusion tensor can only depend on 
the Lorentz force which is quadratic in B, so Kij must be invariant 
under a change of sign of B (Kitchatinov et al. 1994). Thus, we 
have only Kt and kbb that are functions of |B|/i3cq, Pe, and Sc, 
where 



-Boq = ifJ-opU 



2x1/2 
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is the equipartition field strength, and /lo is the vacuum permeabil- 
ity. 

By adding rotation or by imposing a magnetic field, the flow 
becomes anisotropic. The degree of anisotropy is characterized by 
the parameter 



(8) 



2 PRELIMINARY CONSIDERATIONS 

We make use of the fact that in the absence of gravity, and with an 
isothermal equation of state the evolution of entropy is passive, i.e. 
there is no feedback on the momentum equation. We can therefore 
consider the evolution equation for a passive scalar concentration 
per unit volume C, which obeys 



^ = -V ■ (C7C) + ^^V^C, 
ot 



(1) 



where k is the diffusivity and U is the velocity. We note that in the 
compressible case there is an extra term — kV ■ (CV Inp) on the 
right hand side of equation (TJ where p is the density. 

The evolution of the mean concentration C is obtained by av- 
eraging equation lUll, which yields 



dC 
dt 



^ = -V ■ {U C + uc) + kV^C, 



(2) 



where uc = ^ is the mean flux of the concentration density. Under 
the assumption that the spatial and temporal scales of the mean field 
are well separated from those of the turbulence, one can make a 
generalized Fickian diffusion ansatz of the form 



(3) 



where the turbulent diffusion tensor. In the absence of good 

scale separation the multiplication with a diffusion tensor must be 



where and My are the velocity dispersions perpendicular and 

parallel to the direction of anisotropy, f2 or B. In practice, we take 
these directions to be the z direction, in which case u?, — ui and 



is calculated as i (ui 



. Isotropy implies £ = 0. 



3 METHOD 

We simulate turbulence by solving the compressible hydromagnetic 
equations with an imposed random forcing term and an isothermal 
equation of state, so that the pressure p is related to p via p = pCg , 
where Cs is the isothermal sound speed. We consider a periodic 
cartesian domain of size L"^. In the presence of an imposed field 
Bo, we write the total field as i? = Bo + V x A, where A is the 
magnetic vector potential. In the nonmagnetic case, we put Bo — 
0, and with A — initially, we have B = for all times, and the 
hydromagnetic equations reduce to the Navier-Stokes equations for 
p and U, 

°^"^ = -V.t/, (9) 



-2nxU - ctV \np + f + Fvisc + -FLor, (10) 



Dt 
Dt 

where f Lor = in the non-magnetic case, i.e. there is no Lorentz 
force, D/Dt — + [7 ■ V is the advective derivative, fl is 
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the background angular velocity, Fvisc = V ■ 2pvS is the vis- 
cous force, u is the kinematic viscosity, Sij = ^{Uij + Uj^i) — 
^Sij'V ■ U is the traceless rate of strain tensor, and / is a random 
forcing function consisting of plane transversal waves with random 
wavevectors k such that |fe| lies in a band around a given forcing 
wavenumber k{. The vector k changes randomly from one timestep 
to the next. This method is described for example in Haugen et al. 
(2004). The forcing amplitude is chosen such that the Mach number 
Ma = Mrms/cs is about 0.1. 

In the magnetic case, we have a finite Lorentz force f Lor = 
J X B I p, where J = V x i?//io is the current density, /^o is the 
vacuum permeability, and we solve the uncurled induction equation 
in the form 

^ = ?7 X B + ryV'^A, (11) 
at 

where r\ is the magnetic diffusivity, which is assumed constant. 

The test scalar equation is obtained by subtracting equation l[2j 
from equation l[T), which yields 

^ = -V ■ (C7c + uC + uc-uc) +kV^c, (12) 
at 

where c = C ~ C and u — U ~ U are the fluctuating values of 
concentration and velocity and uc is the mean flux of concentration 
density. We restrict ourselves to the case of two-dimensional mean 
fields, C = C{x, z,t), that are obtained by averaging over the y di- 
rection. In the spirit of the test field method, we solve equation il2i 
for a predetermined set of four different mean fields, 

C*^'^ = Co cos fca;, C"" = Co sin fca;, (13) 

Oc Os 

C =Cocosfc2, C —Co sin kz, (14) 

where Co is a normalization factor, that will drop out in the deter- 
mination of Kij since equation ([3) is linear. For each of these mean 
fields C'"' we obtain a separate evolution equation for c^'', 

^ = -V ■ {Ut^" + uC"" + uc"" - HdP^) + kV^c'"', (15) 

where p,q = c, s, or 0. In this way , we calculate four different 
fluxes, = ucPi, and compute the 6 relevant components of 

5 Q 

Ki\ = —{coskxTi —sinkxTi )/k, (16) 

ACis — —{coskzT'i'^ ~ AnkzT'i'^) /k, (17) 

for i = 1, 3. Here, angular brackets denote volume averages. 
The expressions for nn and in equations il6\ and J17I ) differ in 
that the superscripts sO and cO are replaced by Os and Oc, respec- 
tively. Note that the components of Hij for j = 2 are irrelevant for 
mean fields that are independent of y. 

The method provides the functions Kt, ^n, and Knti (or kb 
and kbb) as functions of z and t. Unless otherwise specified, we 
quote averages over z and t. Error bars are calculated by aver- 
aging the results first over z and then over subsequent time in- 
tervals At whose length corresponds to one correlation time, i.e. 
At itrmsfcf ~ 1. This results in a number of independent measure- 
ments from which we calculate the standard deviation, which is 
what we quote as the error. 

The simulations have been carried out using the PENCIL 
COD^I which is a high-order finite-difference code (sixth order 

^ |http : / /Tinni . nordita ■ org/sof tware/pencil-code| 



in space and third order in time) for solving the compressible hy- 
dromagnetic equations. As of revision rl0245, the test scalar 
equations are implemented into the public-domain code similar 
to the way described above and are invoked by compiling with 
TESTSCALAR=testscalar 0. The numerical resolution used 
in the simulations depends on the Peclet number and reaches 256"^ 
meshpoints for our runs with Pe ~ 200. 



4 RESULTS 

As in earlier work on determining the turbulent magnetic diffusivity 
using the test field method, we present our results for the relevant 
components of Kij by normalizing them with respect to the refer- 
ence value 

KtO = ^Urmsk^^. (18) 

This is motivated by the first order smoothing result Kt — ^tu^^^^^ 
and the assumption that the Strouhal number St — rUrmsfcf is about 
unity (Brandenburg et al. 2004). We focus on the case kt/ki ~ 3, 
which is a compromise of having enough wavenumbers between kf 
and the dissipation wavenumber (which is where an inertial range 
can develop), and leaving some room for scales larger than that of 
the energy-carrying eddies corresponding to wavenumber kf. On a 
few occasions we also consider other values of kf/k\. 

We begin by considering the case with rotation and no mag- 
netic field. In Fig. [T] we compare the dependence of Kt, kq, 
and KQQ on Co for fixed values of Pe — 1.8 and 120, using 
Sc = 1 in all cases. Note that Kt shows a mild decline with Co 
while both kq and Knn are positive and increase with Co. For 
Pe — 1.8 the results can be fitted to expressions of the form 
Kt « 0.7/[l + (0.3Co)2], Kn « 0.25Co/[l + (0.5 Co)^], and 
Knn « 0.05 Co/[l + (0.05 Co)^]. For larger values of Pe the de- 
pendence of Kt and Kn on Co is qualitatively similar, but the de- 
pendence of KQ(i on Co now shows indications of a sign change for 
larger values of Co. A positive Kn would agree with the result of 
Riidiger et al. (2005), although their result was explained to be due 
to the additional effect of stratification. 

The dependence of Kt, Kn, and Knn on Pe for Co = 0.5 
and Sc = 1 is shown in Fig. [2] As expected, for very small val- 
ues of Pe all three functions go to zero like Kt/Kto ~ 0.4 Pe and 
Kfi/fCto ~ 0.05 Pe^, while Knn is close to zero and, within error 
bars, essentially compatible with zero. For large values of Pe, both 
Kt and Kn seem to reach asymptotic values, while for Knn there is 
no clear trend. An asymptotic value of Kt was expected based on 
the analogy with the magnetic case where the magnetic diffusiv- 
ity was also found to reach an asymptotic value for large Reynolds 
numbers (Sur et al. 2008). 

In all our cases we find that Kn is positive. This result is in 
agreement with that of Riidiger et al. (2005), who define a quantity 
X such that xQ, = — Kn. They point out that in spherical coordi- 
nates (r, 9, (f)) a finite Kn corresponds to a finite Kr4, which, in turn, 
leads to a finite azimuthal heat flux. Its astrophysical importance is 
limited, because the azimuthal average of its divergence vanishes. 
However, Riidiger et al. (2005) stress that the sign and magnitude 
of Kn provides an important test of mean-field theory. Their re- 
sult was motivated using a result from the mixing-length model 
of Kitchatinov et al. (1994) that Kij should be proportional to the 
Reynolds-stress tensor UiUj. A positive Kn would correspond to a 

^ |http : / /pencil- code ■ googlecode . com/| 
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Figure 1. Co dependence of Kt, ftn. ^nd Kfjn for Pe = 1.8 and 120. 
For Pe = 1.8 the results can be fitted to expressions described in the text, 
while for larger values of Pe the dependences tend to be more complicated. 
The dashed lines indicate the behavior expected from the simple treatment 
of the T approximation in Appendix IaI The slope of the dashed lines in the 
middle and last panel is therefore 1 . For these runs the anisotropy parameter 
e varies between and 0.04. 

negative k^^, which could be explained by a negative UrU^ that 
is expected for stratified turbulence with dominant vertical turbu- 
lent velocities. This argument, which is originally due to Gough 
(1978), does not apply in our case. However, a positive value of 
is an immediate consequence of the Coriolis force producing sys- 
tematic sideway motions; see Fig. [3] That is, suppose there exists 
a mean concentration gradient in, say, the positive x-direction, i.e., 
VC ~ X. Then a small fluid parcel that is moving in the down- 
gradient direction, u ~ — VC, would experience, via the Coriolis 
acceleration, a deflection into the positive j/-direction. The mean 
effect of such deflections would be to produce a positive coher- 
ent flux in the direction given by — VC x Cl, and thus «:n must 
be positive. A positive value of Kn can also be explained using 
the pressureless r approximation which predicts Kn/nt — Co; see 
Appendix |A] This underlines the potential usefulness of the r ap- 
proximation in mean-field theory, even if this approach is used in 
its most rudimentary form. 

The simple treatment of the r approximation in Appendix IaI 
predicts that Kt declines with increasing values of Co, while both 
KQ and Kfin increase and are positive. This is in agreement with 
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Figure 2. Pe dependence for Co = 0.5 keeping Sc = 1 in all cases. The 
resolution is increased with increasing values of Pe and reaches 256"' for 
the run with the largest value of Pe. 




Figure 3. Sketch illustrating the systematic sideway motion caused by the 
Coriolis force and its effect on the flux and hence the sign of kq. 
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Figure 4. Dependence on the magnetic field .strength for Pe = 1.8 (dashed) 
and for Pe in the range 12-15 (solid). For these runs the anisotropy param- 
eter £ varies between and 0.04. 



Figure 5. Pe dependence for Co = and Bq / Bcq = 1.3. The dotted lines 
in the first and last panel correspond to 0.4 Pe and 0.05 Pe'^, respectively. 
For the run with the largest value of Pe we have e = 0.20. 



the simulations, but quantitative details are not reproduced: in the 
simulations the decline of Kt as well as the increase of kq and hqq 
are slower than what is predicted by the simple treatment of the r 
approximation. There could be several reasons for this departure 
of the simulations from theory. Firstly, the definition of Co in not 
unique. However, even for small values of Pe the functions Kt, kq, 
and Knn would require different scaling factors on Co to fit the re- 
sults (see above). Secondly, the feedback of rotational effects on the 
Reynolds stresses are not taken into account in the simple form of 
the r approximation. On the other hand, the rotational anisotropy 
parameter e remains fairly small (below 0.04) even for the largest 
values of Co considered. The correlation fc was neglected, but this 
should be permissible in the current case with (5-correlated forc- 
ing; see Sur et al. (2007) for the corresponding case with magnetic 
fields. Finally, the neglect of the pressure term in the simple form of 
the T approximation is not justified, and a more rigorous treatment 
might alter the results significantly. 

Next, we consider the case with a magnetic field and look at 
the dependence of Kt and kbb on the magnetic field strength for 
fixed Pe=1.8; see Fig.|4] The function kb is also plotted, but, as 
expected, it is compatible with zero. Again, Kt decreases with in- 
creasing field strength, but kbb is finite and positive, and increases 
with Bo- 



The dependence of Kt and kbb on Pe is shown in Fig.|5] Both 
functions increase with Pe as long as Pe is below about 10, and then 
seem to level off, although it would be desirable to confirm this for 
larger values of Pe. The Pe dependence of Kt is similar to that for 
Bo — and Co — 0.6 where Kt goes to zero with decreasing 
Pe like 0.4 Pe, and kbb ~ 0.05 Pe^, which is reminiscent of the 
behavior of Kn for Bo = and Co — 0.6. 

Finally, we consider the dependence of the components of k^j 
on the wavenumber k of the test scalar in equations l |16| l and l |17| l. 
A dependence of Kij on k reflects the fact that there is poor scale 
separation, i.e. k/k{ is not very small. In such a case (as already 
discussed in the multiplication with a turbulent diffusivity in 
equation l|3j must be replaced by a convolution with an integral 
kernel, as in Brandenburg et al. (2008b) for the case of magnetic 
diffusion and a effect. In Fourier space the convolution corresponds 
to a multiplication. As long as we consider monochromatic mean 
fields (proportional to e''° "' for a single wavevector) the values of 
the numerically determined components of k^j correspond to the 
value of the Fourier transform kij at wavevector k. The full integral 
kernel can then be assembled by determining the full k dependence 
and then Fourier transforming back into real space. 

In Fig. |6]we show the resulting dependences Kt(fc), Kn(fc), 
and Knn (k), similar to earlier work on turbulent magnetic diffusion 
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Figure 6. Wavenumber dependence of Kt , kq. and fcnn for Co = 0.6 
and fcf/fci = 8 with Pe = 70 compared with the case fef/fci =3 with 
Pe = 17. In the upper two panels the lines represent curves of the form 
119) with a = 0.62 and 0.5, respectively. For these runs the anisotropy 
parameter e is around 0.003 for the runs with fcf/fci =8 and Pe = 70 and 
around 0.036 for the runs with k{/k\ =3 and Pe = 17. 



(Brandenburg et al. 2008b). It turns out that both Kt and kn can be 
fitted reasonably well to a Lorentzian, 

" 1 + {ak/hY ' ^^'^^ 

where a ~ 0.62. This value of a is close to the corresponding 
value for turbulent magnetic diffusion where a ~ 0.5 was found 
for the isotropic case (Brandenburg et al. 2008b) and a ~ 0.7 in 
the anisotropic case with shear, but no rotation (Mitra et al. 2009). 
However, the coefficients Kto and Una depend on value of kf/k\. 
For fcf/fei = 3 we find Kto/'^to = 1-4 and kno ~ 0.3, while for 
fcf/fci = 8 we find Kto/'i^to = 1-5 and ftno/^^to ~ 0.5. To make 
the curves for different values of k{ and Pe more nearly overlap, 
we have scaled Kt by a factor St = kt/ki, kq by a factor sn = 
(fef/fci)^'^, and Knn by a factor sn = (kf/ki)"'^. 

A surprising result occurs for knn (k) which takes large neg- 
ative values for k/k{ ~ 1, but tends to vanish for k 0. This 
means that the diffusion of mean fields along the axis of rotation is 
strongly suppressed once its typical length scale becomes compa- 
rable with the scale of the energy-carrying eddies. 



Figure 7. Dependence of Kt, kq, and Knn on the forcing wavenumber 
fcf/fci for A: = with Co = 0.6 and Pe ks 17. Note the scaling of Kt and 
Kn with fc r ^ , while Knn is essentially compatible with zero. 



The two cases shown in Fig.|6]differ not only in the value of 
fcf / ki (3 and 8), but also in the value of Pe (17 and 70). In order to 
show that the dependence of the scaling coefficients Kto and Kno is 
primarily due to the change in fcf/fci we plot in Fig.|7]the depen- 
dence of Kt, Kn, and Knn on the forcing wavenumber fcf/fci for 
fc = with Co — 0.6 and Pe ^ 17. It turns out that both Kt and 
Kn decrease with fcf like fcj^^, while Knn is essentially compatible 
with zero. 

In order to illustrate the meaning of the error bars, we show 
in Fig. [8] time series of Kt, Kn, and Knn for a run in Fig. |6] with 
fc = 0, fcf/fci = 8, Co = 0.6, and Pe — 70. Time is given in 
turnover times and covers nearly 200 in those units. This results in 
nearly 200 independent measurements, whose standard deviation 
yields the one-sigma probability range of finding measurements to 
lie in that interval. 

In Fig. |9] we show the k dependence of the transport coeffi- 
cients in the presence of a magnetic field and without rotation. The 
k dependence of Kt is about that found in Fig. |6l but now kbb 
shows a marked dependence with a sign change at k/kf « 0.5 
about values of kbb of opposite sign with extrema of similar mod- 
ulus, KBB/wto ~ +0.3 at fc/fcf !^ and kbb/hw ~ —0.3 at 
k/kf ~ 1. Furthermore, for fc ^ and fc ^ oo the values of kbb 
are small. 
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Figure 8. Time series of Kt, kq, and kqh for a run in Fig.|6]witli k = 0, 
fcf/fci = 8, Co = 0.6, and Pe = 70. The average values are shown as 
straight soUd lines and the en'or margin is indicated by dashed lines. 



Figure 9. Wavenumber dependence of Kt, kb, and kbb f™ Co 
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5 CONCLUSIONS 

The test scalar formalism described above proves to be a useful 
tool for calculating the components of the turbulent diffusion ten- 
sor for the mean passive scalar concentration. Here the mean pas- 
sive scalar concentration is defined as a two-dimensional field ob- 
tained by averaging over one coordinate direction. Therefore the 
diffusion tensor has 6 relevant components. In the present inves- 
tigation, we have only looked at the case of homogeneous turbu- 
lence where anisotropy is produced either by rotation or by a mag- 
netic field. Another interesting alternative would be to look at the 
effects of a linear shear flows of the form = (0, Sx, 0). In- 
vestigations of the corresponding magnetic case with the test field 
method have already been carried out (Brandenburg et al. 2008a), 
and they confirmed the similarity between cases with shear and ro- 
tation. The application of the test scalar method to cases with shear 
would therefore be straightforward. 

The test scalar method is analogous to the test field method 
for obtaining turbulent transport coefficients for the evolution of 
the mean magnetic field. Therefore, many of the questions raised 
in that case carry over to the present case. An example is the de- 
pendence of the transport coefficients on the length and time scales 
of the mean field. Usually, rapid spatial and temporal variations are 
less efficiently diffused by turbulence whose correlation length and 



time are no longer small compared with the corresponding scale of 
variation of the mean field. 

A surprising result of the present work is the fact that turbulent 
transport along the rotation axis is not strongly reduced unless the 
length scale of the mean field becomes comparable with the corre- 
lation length of the turbulence, i.e. when A; ~ fcf . In that case, there 
is significant suppression of the turbulence along the direction of 
rotation. A similar result also applies if anisotropy is produced by 
an applied magnetic field. 

The present investigation is restricted to the case when the 
scalar field is passive. This applies to the evolution of the entropy 
density in situations where an isothermal equation of state applies 
and where gravitational stratification is negligible. However, once 
gravitational stratification or non-isothermal effects become impor- 
tant, the test scalar formalism can only be used to describe small 
departures from a given reference state. It would be worthwhile ex- 
ploring this case further 

Our work has shown that for non-helically forced isotropic 
turbulence, i.e. without rotation or magnetic fields, the turbulent 
diffusivity of a passive scalar is similar to the value of the turbu- 
lent diffusivity of the magnetic field in the kinematic limit. While 
this was expected based on results from the first-order smoothing 
approximation, the situation is less clear for turbulent kinematic 



© 2009 RAS, MNRAS 395. [T399iri606l 



1606 A. Brandenburg et al. 



viscosity, which suggest that the turbulent magnetic Prandtl num- 
ber should be 2/5. Simulations of Yousef et al. (2003), based on the 
decay rates of large-scale magnetic and velocity structures, have 
shown that the turbulent magnetic Prandtl number is indeed also 
unity. However, their method is not very accurate and it would be 
good to reconsider this by developing a method analogous to the 
test field and test scalar methods for test flows. 
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APPENDIX A: THE POSITIVE SIGN OF Kn 

In this appendix we motivate the origin of the positive sign of kq 
using the pressureless r approximation. In this approach one con- 
siders the evolution equation for ,F, 



dt 



iiiC + UiC, 



(Al) 



where dots denote time derivatives (Blackman & Field 2003, Bran- 
denburg et al. 2004). We use equations dlOb and l ll2b . assume that 
u — U, i.e. that U — 0, ignore the pressure term and magnetic 
fields, assume /c = 0, and obtain 



iiiC = —2Qj€ijkUkC + triple correlations, 
UiC = —UiUj y jC + triple correlations. 



(A2) 
(A3) 



The triple correlation terms result from the nonlinearities in the 
evolution equations equations l llOt and l ll2t . In the r approximation 
we substitute the sum of the triple correlations by quadratic correla- 
tions, i.e. by —Uic/r in the present case (Vainshtein & Kitchatinov 
1983, Kleeorin et al. 1990). We write the resulting equation in ma- 
trix form. 



Qjr^ _ _ 

T^r^ = -\-ikTk - TUiUj VjC, 

at 



(A4) 



where L^fc = 5ik + 2Q,jTeijk- We solve this equation for ^ and 
obtain 



Ti = -(L )ij TUjUk VfeC + T 



dt 



where 




1 Co ^ 
Co 1 

1 + Co^ 



(A5) 



(A6) 



(A7) 



In the stationary state we may ignore the time derivative and re- 
cover equations lO and Q with 



Co 



2 ' "~ = Co, 



Co^+e 



(A8) 



where e — 1 ~ u^/u^^ was defined in equation JS}. In the cases 
reported here the value of e was never found to be negative. Note 
also that 1 — e — /ti^ in the denominator of equation l lASt is 
always positive. Therefore, in addition to nn being positive, also 
Knn should be positive. 
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